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SUMMARY 

This  report  presents  a  unified  treatment  of  theory  for 
calibrating  the  six-port  ref lec tame ter  and  the  dual  reflectometer 
network  analyser  for  measuring  power  and  scattering  coefficients, 
given  the  minimum  of  standards  for  each.  It  describes  a  reflecto¬ 
meter  in  terms  of  voltage  variables  and  reduces  the  calibration 
theory  to  equations  in  real  unknowns  to  facilitate  encoding  the 
computation  necessary  for  known  calibration  methods  in  a  low-level 
computer  language.  The  treatment  forms  the  basis  for  software  in 
use  with  six-port  ref lectometers  and  dual  six-port  network  analysers 
in  development  at  RSRE  for  comparing  RF  and  microwave  metrology 
standards . 
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INTRODUCTION 


Measurement  of  power  and  complex  voltage  transmission  and  reflection 
coefficients  with  a  six-port  ref lec tome ter  were  first  described  by  Hoer  and 
Engen  more  than  a  decade  ago  (1,2).  This  form  of  ref lectometer  consists 
essentially  of  four  radiation  detectors,  each  indicating  the  power  absorbed  by 
it,  arranged  around  a  structure  guiding  radiation  to  the  detectors  and  to  the 
device  under  test  (DUT) .  The  waveguide  structure  serves  solely  to  direct  to 
the  detectors  different  samples  of  waves  incident  on  and  reflected  from  the  DUT 
and  observations  are  made  of  the  ratios  of  power  absorbed  by  three  of  the 
detectors  to  that  absorbed  by  the  fourth.  This  simplicity  of  principle  implies 
that  this  form  of  ref lectometer  may  be  used  for  measuring  scattering  coefficients 
at  any  frequency  at  which  artefacts  can  provide  physical  approximation  to  a 
conceptual  impedance  for  calibrating  the  instrument.  Examples  have  been 
reported  operating  over  the  10  MHz  to  1  GHz  range  by  Woods  and  Granville-George(3) , 
from  2  to  18  GHz  by  Hoer  and  Engen  (4),  and  up  to  75  GHz  by  Weidman  (5). 

This  detailed  report  presents  a  unified  treatment  of  the  theory  of  known 
procedures  for  calibrating  both  the  six-port  ref lectometer  and  the  dual  reflecto- 
meter  network  analyser  (6-8).  Because  on-line  computation  of  results  is  an 
operational  necessity,  the  theory  has  been  reduced  to  equations  in  real  unknowns 
without  ambiguity  of  square-root  signs  to  facilitate  programming  in  a  low-level 
computer  language.  The  treatment  was  developed  as  a  necessary  prerequisite  for 
programming  desktop  computers  controlling  these  instruments  for  the  eventual  aim 
of  comparing  standards  for  RF  and  microwave  metrology.  To  meet  this  aim  the 
calibration  methods  treated  are  mainly,  but  not  exclusively,  those  requiring  the 
minimum  of  known  standards. 

Although  the  operation  of  six-port  instruments  depends  on  the  observation 
of  power  ratios,  we  start  from  a  description  of  an  n-port  ref lectometer  in  terms 
of  voltage  variables.  This  is  appropriate  because  a  scattering  coefficient  is 
merely  a  ratio  of  two  complex  numbers  each  representing  a  voltage.  This 
description  also  has  the  merit  of  enabling  the  relatively  few  assumptions  under¬ 
lying  six-port  ref lectometer  use  to  be  identified.  Such  identification  is  as 
necessary  for  quantifying  uncertainty  of  measurement  with  a  six-port  as  it  is 
for  any  other  form  of  ref lectometer. 

Starting  from  a  description  of  the  n-port  reflectometer  also  facilitates 
establishing  equivalence  relations  between  six-  and  four-port  ref lectometers. 

These  relations  enable  two  six-port  ref lectometers  forming  a  network  analyser  to 
be  calibrated  in  terms  of  a  reflectionless  waveguide  plus  a  reflecting  termina¬ 
tion  of  unknown  but  repeatable  voltage  reflection  coefficient  (VRC).  This 
method  (8)  calibrates  both  ref lectometers  in  terms  of  the  impedance  ZQ  character¬ 
ising  the  mode  transmitted  by  a  known  length  of  uniform  waveguide  (reflectionless 
by  definition)  and  reflected  by  the  termination.  It  yields  the  propagation 
constant  of  the  waveguide  and  if  this  be  air-dielectric  transmission  line  then 
subsequently  measured  results  can  be  normalised  to  the  nominal  real  RQ  (af^Z0) 
characterising  the  TEM  mode  postulated  for  a  lossless  transmission  line. 

2  INTRODUCTION  TO  SIX-PORT  REFLECTOMETER  CALIBRATION 

2.1  THE  GENERAL  N-PORT  (n  >  3)  REFLECTOMETER 

We  define  a  reflectometer  as  an  n-port  linear,  but  not  necessarily 
reciprocal,  time- stationary  waveguide  junction  (where  n  >  3)  directing 
radiation  from  a  source  at  port  l.  Figure  1,  to  an  output  (measurement) 


FIG.  1.  GENERAL  n-PORT  REFLECTOMETER  (n>3 ) 


port  2  and  to  detectors  at  ports  3  to  n.  We  assume  that  there  is  one  mode 
transmitted  at  each  port  but  not  necessarily  the  same  mode  at  all  ports  and 
that  there  is  no  evanescent  mode  at  any  port;  these  assumptions  enable  us 
to  include  mode  transformers  within  the  junction.  They  also  enable  us  to 
define  at  the  reference  plane  of  each  port  i  (where  n  >  i  >  0)  voltages 
proportional  to  the  electric  field  intensity  described  by  the  transmitted 
mode  of  waves  incident  on  and  emergent  from  port  i.  We  denote  these  by 
complex  numbers  aj,  and  b£,  respectively ,  such  that  their  RMS  amplitudes  are 
|a£ 1  and  |b^|.  The  assumption  of  junction  linearity  implies  that  each  b£ 
is  composed  of  simply  the  sum  of  contributions  each  proportional  to  an  a£ 
and  this  can  be  conveniently  represented  by  the  matrix  notation 

[b]  -  [S][a]  (2.1) 

in  which  the  elements  of  each  row  of  l  S]  are  complex  numbers  relating  the 
corresponding  b£  of  the  column  matrix  { b]  to  every  a£  forming  the  column 
matrix  [a].  The  elements  of  IS]  can  be  considered  as  scattering  coefficients 
normalised  to  a  yet-to-be-def ined  impedance. 

Without  defining  this  impedance  we  can  define  for  each  detector 
connected  to  port  i  (where  now  n  >  i  >  2)  a  VRC  r*  -  a£/b£.  These  can 

be  substituted  in  equation  (2.1)  to  give 

*1  “ 
b2 
b3 

b4  -  IS] 

i 
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For  n  >  i  >  2  this  notation  represents 


..  0  - 
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We  now  make  the  further  assumption  that  there  is  unequal  transmission 
between  the  ports  of  the  junction  (a  condition  hard  to  avoid  in  practice). 
This  implies  that,  denoting  the  square  matrix  of  equation  (2.2)  by  [fil ,  the 
elements  of  any  row  of  (Hi  are  not  in  general  multiples  of  the  corresponding 
elements  of  any  other  row  so  that  dat[n)  *  0.  Thus  [H)  may  in  general  be 
inverted  to  give  [Hi  “  [ Hi  and  we  can  write 


-  [HI 


0 

0 


or  its  equivalent 


*1  ■  nUbl  *  m12b2 
*2  *  "2lbl  *  “22b2 
bi  *  "ilbl  +  mi2b2  <for  “  *  1  >  2) 


(2.3.1) 

(2.3.2) 

(2.3.3) 


where  m££  and  m£2  denote  elements  of  the  first  and  second  columns  of  [HI, 
respectively. 

Rewriting  equation  (2.3.2)  as  bj  ■  (a2-m22b2)/m21  and  substituting  in 
equation  (2.3.3)  gives 


Defining  T  =  a2/b|  *8  tlie  VRC  of  an  artefact  connected  to  port  2  (either  a 
DOT  or  a  calibration  standard)  allows  us  to  write  from  this  equation 


or 

b.  -  («i  r  ♦  e.)b2  (2.5) 

where  ■  m^/n^i  end  &£  -  (m£2m21  "  ®22®il) /m2l  are  dimension-less 
constants  characterising  the  waveguide  junction. 

Equations  (2.4)  and  (2.5)  are  fundamental  to  all  ref lectometers 
comprising  a  time-invariant  waveguide  junction.  Their  application  to 
varying  junctions  is  dependent  on  the  validity  of  additional  assumptions 
such  as  the  repeatability  of  RF  switches  within  the  junction.  Application 
of  either  equation  to  a  particular  ref lectometer  depends  on  the  method  of 
detection  employed  since  it  is  the  detector  output,  or  some  function  of  it, 
that  forms  the  observed  indication.  In  the  next  section  we  apply  these 
equations  to  a  four-port  ref lectometer  to  illustrate  this  dependence  on 
detector  output  and  to  establish  an  equivalence  between  the  four-port  and 
a  perfect  ref lectometer  preceded  by  a  two-port  "error  box",  a  representation 
we  later  use. 

2.2  THE  FOUR-PORT  REFLECTOMETER 

For  this  n  «  4  in  Figure  1,  and  the  two  detectors  D3  and  D4  are  each 
assumed  to  provide  an  output  voltage  V£  (i  ■  3,4)  proportional  to  the  RF 
voltage  presented  to  it  by  the  four-port  waveguide  junction.  These 
voltages  Vj.  are  at  an  IF  to  enable  both  the  real  and  imaginary  components 
of  v^/v^  to  be  indicated  by  an  associated  superheterodyne  receiver  (9). 

We  initially  assume  that  the  complex  indication  provided  by  the  receiver  is 
some  function  of  this  ratio,  ie  indication  w  ■  fn(v4/v3).  We  relate  w  to 
the  b3  and  b4  of  equation  (2.5)  by  first  writing 

vi  "  **i+bi>  *  (i+ri)bi 

and  substituting  in  (2.5)  to  give 

vi  "  <1+ri>(air+6i>b2 

V  -  fn(v4/v3)  -  fn((l+r4)(a4r+eA)/((l+r3)(a3r+e3)))  (2.6) 

or  w  »  fn((dl>e)/(cr+l) 

where  c  »  <*3/83  »  d  •  (l+r4)o4/((l+r3)63) 

e  -  u*r4)34/((ier3)  e3) 

ere  characteristic  of  the  ref lectometer. 


5 


If  the  additional  assumption  be  made  that  the  receiver  indication  w  is 
a  linear  function  of  V4/V3,  absorbing  the  coefficients  of  linearity  in  c,d,e 
equation  (2.6)  becomes  the  bilinear  transformation  between  w  and  T  given  by 

w  -  (dr+e)/(cr+l)  (2.7) 

Physical  significance  is  imparted  to  w,  enabling  T  to  be  found  from  equation 
(2.7),  only  by  calibrating  the  instrument.  Because  it  is  characterised  by 
three  complex  numbers  c,d,e,  it  can  be  calibrated  by  observing  vj,(f  ■  1,2,3) 
corresponding  to  the  connection  of  three  standards  providing  different  known 
r^.  Therefore  it  is  the  impedance  in  terms  of  which  these  1*4  are  known  that 
defines  the  Z0  to  which 

i.  the  result  of  measurement  T, 

ii.  the  elements  of  [  S]  in  equation  (2.1),  and 

iii.  the  of  the  detectors 

are  all  normalised.  It  is  this  ZQ  that  determines  c,d,e  in  the  bilinear 
transformation  (2.7)  to  completely  describe  the  ref lectometer  (comprising 
the  four-port  junction,  detectors  and,  given  linearity  of  response,  the 
associated  receiver). 

The  explanation  in  the  preceding  paragraph, that  the  impedance  left 
undefined  in  section  2.1  is  defined  by  the  ZD  to  which  the  r&  of  known 
calibration  standards  is  normalised,  applies  equally  to  any  ref lectometer. 

It  will  be  shown  in  section  3  that  a  modification  of  the  calibration 
method  just  outlined  can  be  applied  to  the  six-port  ref lectometer,  so  we 
proceed  to  derive  theory  appropriate  to  four-port  ref lectometer  calibration 
for  later  use.  Writing  for  at  least  three  known  VRC  r*  *  x4  +  jya 
(j  =  /-T;  l  ■  1,2,3)  and,  for  the  corresponding  wjj,  ■  +  jv^,  and  denoting 

the  real  and  imaginary  components  of  c,d,e  by  suffixes  R  and  I  respectively, 
then  equation  (2.7)  becomes 

<u*  +  jv4){l+(cR+jCI)(*i+jyf>}  "  (dR+jdI)(xt+jyA)  +  (eR+^*I)  (2'8) 

Separating  real  from  imaginary  components  gives  for  i  ■  3  an  exact,  and  for 
i  >  3  an  overdetermined,  set  of  simultaneous  linear  equations  in  the  six 
real  numbers  constituting  c,d,e: 


u* 

v* 


dRXrdiy£  +  eR-  (**U£-yiV  CR+  <Vi  +  y*U*)cl) 
dRy£ +  dixz  +  ei~  ^x£v£ +  ~  ■  yjtV0! ) 


(2.9) 


After  equations (2.9)  have  been  solved  for  c,d,e  (using  a  least-squares  method 
if  more  than  the  minimum  of  standards  are  used) ,  they  can  be  recast  to  allow 
the  VRC  of  a  subsequently  connected  DUT  to  be  found  from 


x  -  (p(u-eR)  -  q(v-ej)}/(p2+q2) 
y  -  (q(u-«fl)  ♦  p(v-Bj) )/ (p2+q2) 


(2.10) 


where  p  -  (dR-ucR+vCj)  and  q  -  (-dj+vc^uCj) 
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Thus  any  unknown  r  ■  x+jy  can  in  principle  be  found  from  the  complex 
indication  v  provided  that  p  and  q  are  not  both  zero  (which  we  show  later 
to  be  true  generally).  The  words  "in  principle"  have  been  used  because 
measurement  relies  on  (2.7)  being  used  to  interpolate  in  both  x  and  y 
directions  in  the  complex  r  plane  between  the  known  r*  ■  x^+jyj  provided 
by  the  standards;  extrapolation  beyond  those  values  may  not  be  valid 
(because,  for  example,  the  detector  response  may  no  longer  be  linear). 

An  alternative  treatment  of  the  four-port  reflectometer  regards  it  as 
equivalent  to  an  "error  box"  (ie  a  two-port  junction  having  equivalence  to 
c,d,e),  one  port  of  which  forms  the  measurement  port  and  having  connected 
to  the  other  a  perfect  reflectometer  (one  whose  indication  w  is  simply  a 
VKC)  as  shown  in  Figure  2. 


0 

Source 


Notional 

perfect 

reneciomeier 


Measurement 

port 


FIG.  2.  ERROR  BOX  REPRESENTATION  OF  REFLECTOMETER 


Ve  derive  this  equivalence  for  later  use  by  remembering  that  a  two-port 
junction  with  scattering  matrix 


transforms  a  VKC  T  into  w,  where 


w 


’ll 


.  *12*21 


r 

r 


*12*2l\  +  *12*21 

*22  /  *22^“*22r) 
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Comparing  this  with  aquation  (2.7)  ravrittan  aa 


e» 


f 


V 


d  .  e-d/c 
c  cT  ♦  1 


(2.1i) 


shows  tha  equivalence  of  tha  "error  box"  raprasantation  to  ba  valid  if 

e  -  -s22  ,  d  -  *12*21  “  *11*22  *  *  "  *11 

Both  reprasantations  account  for  RF  signal  leakage  between  the  detectors, 
since  transmission  between  all  ports  of  tha  waveguide  junction  was  initially 
postulated;  neither  takes  account  of  any  interdependence  of  phase  and 
amplitude  within  the  associated  receiver.  (We  remark  that  they  may  be 
applied  to  tha  reflec tome tars  of  a  dual  ref lectometer  network  analyser  (see 
section  4)  only  if  there  are  no  leakage  paths  between  the  ref lactometers.) 


The  error  box  representation  makes  it  evident  that  the  ref lectometer 
cannot  function  if  si2®21  "  0*  We  ®*y  therefore  assert  that  812*21  *  ®  aiu* 
this  is  equivalent  to  d-ce  *  0.  But  from  equation  (2.10) 

P  -  jq  -  dR  +  jdj  -  (Cn+jc!) (u+jv) 
so  that  p2  +  q2  •  !p-jq|2  ■  |  d-wc  | 2 

Hence  the  condition  p  «  q  ■  0  is  equivalent  to  w  *  d/c  in  (2.11)  which  is 
inadmissible  since  d-ce  *  0.  Thus,  as  was  asserted  earlier,  T  “  x+jy  can 
always  be  found  from  w. 

2.3  THE  SIX-PORT  REFLECTOMETER 


For  this,  n  »  6  in  Figure  1,  and  each  of  the  four  detectors  D3  to  Dg 
is  assumed  to  provide  a  linear  indication  of  the  power  it  absorbs.  (In 
practice  bolometric  detectors  are  most  frequently  used  because  of  their 
excellent  approximation  to  this  assumption  (10)).  If  each  indicates  power 
?i  -  Gi(|bi]*-  |a£ | 2)  -  Gi(l-|r$J2)  |b£ | 2,  (where  6  >  i  >  2)  then  Gj  may  be 
interpreted  as  tha  conductive  component  of  an  admittance  characterising 
the  mode  transmitted  at  port  i.  From  equation  (2.4)  we  form 


\\\‘ 


so  that 


Pi  " 


VHrJ  > 


L  vu' 

V12  “21  I 

1  i**i* 

r  ♦  L 

m22mil\| 

Ki  \i2 

a 

N> 

1-^ 

b,r 


♦  e 


ib2r 


where  <»£  and  8£  are  defined  by  0£2  •  G^d—  | T£ | (mn/m2i)‘ 

and  S£2  -  Ci(l-|ri|2H(mi2m21  -  ®22mil)/m2l)2 

by  analogy  with  the  definitions  at  (2.S). 


(2.12) 
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r 


Defining  |wjJ  •  P^/P^  for  6  >  i  >  3  end  k  ■  i  -  3  gives 

12 

l»kr 


or 


l»J' 


(2.13) 


where 


°3/63 


°l/63  • 


V63 


Fran  (2.13)  the  three  observed  power  ratios  can  be  interpreted  as 
defining  three  circles  in  the  complex  plane  intersecting  in  r  ■  x+jy  and 
operation  of  the  six-port  ref lectometer  depends  on  this  intersection  being 
unique,  for  which  the  centres  must  be  neither  coincident  nor  collinear. 

The  centres  can  be  found,  for  dropping  the  suffix  k,  equation  (2.13)  can 
be  written  as 


|w|2(cr+l)(c*r*+l)  -  (dr+e) (d*r*+e*) 

where  *  denotes  complex  conjugate. 

|w|2(|c|2|r|2+cr+c*r*+i)  -  |d|2|r|2+de*r+d*er*+|e 


(2.14) 


| r 1 2 _ c- r  .  -,S*M2-d*e  r 

I*  I  .  .1  .  1  .  .  .9.  .*>  1 


I  2  ' 

|w|  -|e 


nr-i 


C  w 


m -icrw 


|d|2-|c|2|w|2 


9  A  O  O 

Expressing  this  in  the  form  |r|  -r«*r-T«r*  *  r  -|r0|  ,  where  r  is  positive 
real,  identifies  it  as  the  circle  | r— r0|2  “  r2.  Its  centre  is 


r  - 


c*|.|2-d*. 


(cR|w|2-(dRV'>iei)>^(<diVVi)‘cil“|2) 


°  l<i|2-|d|2|»|2  |d|2-  |c|2|w|2 

Since  |ro|2«roro*,  it  follows  from  the  above  that 


<|d|2-|c|2|v|2)2|rj2  -  (c|w|2-de*)(c*|»|V.) 

*  |d j2|e|2+|c j2|w|^-(d*ce+dc*e*) |w|2 

and 

<|d|2-|c|2|«|2)2(r2-|r0|2)  -  (|w|2-|e|2)(|d|2-|c|2|w|2) 

-  -|d|2|e|2-|c|2|w|*+(|d|2+[c|2|e|2)  |w|2 

Adding  gives 

2  B  Iwl 2(|d | 2+ |c | 2 I e |2-d*ce-dc*e* )  _  lw{ 2 |d  -  ce|2 

(|d|2-|c|2|.|2)2  (|d|2-|c|2|w|2)2 
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Thus  for  each  k,  |«^J  generates  the  circle 
/  ~(dReR'*'dlel)k\  /  \ 

V  UJ2-M2|w.!2  /  v  |dkl2-|c|2|wj2  / 


iv^iAkl! 

(|dkl2-|c|2!wk|2)2 


(2.15) 

Now  the  assumption  of  unequal  transmission  between  all  ports  implies  that 
dk  and  e^  are  different  for  all  k  and  equation  (2.15)  shows  that  the  circle 
centres  are  therefore  neither  coincident  nor  collinear.  Generally  both  the 
centre  and  radius  of  each  circle  is  dependent  on  |w|  but  if  c  *  0  then  the 
radius  of  each  circle  is  simply  r^  ■  | w/ d | ^  and  the  centre  co-ordinates 
(Tq)^  are  independent  of  |w|.  It  follows  from  (2.12)  and  (2.13)  that 

c  *  0  implies  that  detector  D3  is  isolated  from  the  wave  a£  reflected  onto 
the  measurement  port  so  that  equation  (2.13)  becomes  a  linear  relation 
between  |wk|  and  T;  this  has  been  utilised  for  one  form  of  five-port 
ref lectometer  (11). 


In  general  then,  T  ■  x+jy  is  uniquely  determined  from  the  three  circles 
each  dependent  on  the  power  ratios  IwjJ^  and  can  be  found  if  sufficient  can 
first  be  found  of  the  seven  complex  numbers  c,  d^,  e^.  One  scheme  for 
calibrating  a  six-port  for  measuring  VRC  has  been  described  by  Woods  (12). 

It  requires  the  connection  of  seven  standards  providing  different  known 
Tj,  (1  <  i.  <  7)  and  observing  the  twenty-one  corresponding  |w^|j  but  its 
theory  is  particularly  simple  to  derive  from  equation  (2.14).  Noting  that 

cT  ♦  c*r*  •  2(cRx  -  Cj.y) 


and 

dre*  ♦  d*r*e  -  2(dReR+diei)  x  -  2(dIeR-dReI)y 
gives  (2.14)  in  the  form 

| w 1 2  (1*  I  c  I  ■ 2 1  r  1 2+2 (cRx-Cly)  )  -  I  d  1 2 1  r  1 2+ 1  e  1 2+2  (d^+d^j)  x-2  (dJ.eR-dReI)y 
Reintroducing  suffix  k  and  writing 

■V  ■  <dRVdIeI)k  and  \  ■  -(dlVdR*I) 

gives 

I  »k  1 1  a*  I  c  1 2  (X^y J)  ■ *2  (Cr Vc  jyt) )  -  Uk|2(*jj*y2H«kl  Vl,2Vl 


(2.16) 


Observing  |w|Jj[  gives  for  each  k  an  exact  set  of  seven  simultaneous  linear 
equations  (2.16)  in  real  unknowns  to  enable  |cp,  cr,  cj,  |djJ2,  lekl2»  Pk» 
qk  to  be  found.  If  more  than  seven  standards  are  available  then  the  over¬ 
determined  set  of  equation  (2.16)  can  be  solved  by  standard  methods  (see 
section  6).  The  solutions  can  then  be  used  in  equation  (2.16)  recast  as 
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(*2+y2) ( |c  J 2 1 vk| 2- |dk| 2)+2(cR| wk| 2-pk)x-2 (cx |wfc 1 2+qk)y  -  I  \  1 2_ I wk  I 
to  find  T  ■  x+jy  from  the  corresponding  observed  |wjJ2.  Writing 


Pk“CR<Wkl 

|c|2|wk|2‘ldk 


qk  +  cllwkl2 

,2  »  8k  "  |  1 2 1  ,2  ,  .  i2  » 

I c I  |vkl  -|dk| 


ek|2"lwkl2  .2  2 

2,  ,2  k  ,2+fk+  8k 


K  icnwkr-idki 

finally  gives 

(*-fk>2  *  (y-g,,.)2  -  r2 
which  can  be  solved  by  writing 


£  - 

giving 


2  2  .^2 .  2.  ,.2A  2, 
rrr2+(f2+82)‘(fl+8l) 
2(f2^f1) 


(83~'82  ®2**8l\ 


F  - 


2  2  ,,2  2N  te2  2N 
r2"r3  (f  3  83J “ 2  g2) 
2(f3-f2) 


and  x  ■  C 


-esj- 


Other  methods  of  calibration  have  been  described  which  rely  on  the  use 
of  only  four  standards  namely  four  short-  or  open-circuits  (13).  We  proceed 
in  the  next  section  to  derive  a  modification  of  the  method  derived  in  section 
2.2  for  a  four-port  to  enable  a  six-port  ref lectometer  to  be  calibrated  by 
means  of  only  three  precisely-known  and  one  approximately-known  standard. 

This  will  pave  the  way  to  treating  the  dual  six-port  network  analyser  and 
its  calibration  in  terms  of  just  one  known  standard. 

3  SIX-  AND  FOUR-PORT  REFLECTOMETER  EQUIVALENCE 
3.1  SIX-  TO  FOUR-PORT  REDUCTION  THEORY 

The  six-  to  four-port  reduction  algorithm  treated  in  this  section  is 
based  on  a  derivation  given  by  Engen  (6).  It  enables  a  vector  "indication" 
to  be  calculated  from  the  power  ratios  observed  for  a  six-port  equivalent 
to  that  given  by  the  receiver  of  a  four-port  reflectometer. 

We  define,  consistent  with  (2.13),  wk  -  (dkr+ek) /(cT+i)  and  solving 
for  T  gives 

r  -  -  *2~*2  .  *3 "  W3 

cwj-dj  cw2-<,2  cw3_d3 


11 


Writing  Qfc  - 

and 

eliminating  f  produces 

the  three 

equations 

«1  - 

l",!2 

• 

f 

a2q2 

2> 

-  iw^r 

b2q3  - 

jwj-n|2  (3.1) 

with 

A  * 

cVdi 

B  - 

c*i-di 

;  m  ■ 

■'lWi 

dle3~d3e 

ce2-d2 

1 

ce3-d3 

ce2~d2 

’  °  ce3-d3 

These  Chree  equations  in  w^  represent  three  circles  in  the  complex  w4 *  Uj+j\^ 
plane  intersecting  in  wj,.  Their  intersection  is  unique  to  each  T  because 
the  three  centres  (0,mtn)  are  bilinear  transformations  of  those  of  (2.14). 
Now  Ok  are  the  observed  power  ratios  and  since  eliminating  T  has  shown  that 
A  and  B  are  constant  for  all  r,  it  follows  that  they  are  characteristic  of 
the  ref lec tome ter.  They  represent,  in  the  w^  plane,  ratios  of  the  radii  of 
the  two  circles  centred  at  (m,n)  to  that  of  the  circle  centred  at  the 
origin.  We  proceed  to  define  three  more  real  constants  characterising  the 
ref lectometer,  namely  the  distances  between  the  centres  of  these  three 
circles  in  the  complex  wi  plane. 

Writing  m  ■  Mcosu  +  jMsiny,  n  *  Ncosv  ♦  jNsinv,  w4“  u^  +  jv^  (M,N  >  0) 
into  (3.1)  gives 


2  2 
Ql  -  u1  +  v, 

A2Q2  •  (Uj-Mcosy)2  ♦  (Vj-Msiny)2 

B2Q3  «  (Uj-Ncosv) 2  +  (Vj-Nsinv)2 

Now  define 

p  *  |m-n|2  ■  M2  ♦  N2  -  2MNcos(v-y)  ; 
q  3  |n|2  -  N2  ; 
r  =  |m| 2  -  M2 


(3.2) 


(3.3) 


These  definitions  of  m,n  associated  with  equations  (3.1)  are  independent  of 
F  so  that  these  distances  between  the  circle  centres  in  the  w^  plane  (p,q,r) 
also  characterise  the  ref lectometer  (the  p  and  q  so  defined  are  not  related 
to  those  of  (2.10)  or  (2.16)). 


All  five  positive  real  numbers  A^,B^,p,q,r  can  be  found  by  observing  the 
power  ratios  corresponding  to  nine  (or  more)  OUT  whose  VRC  are  known  only  to 
be  different;  their  precise  values  need  not  be  known.  W#  write  from  equations 
(3.2)  and  (3.3) 


A  ■  -2Mu1cosy-2Mv1siny+M*  and  B*Q,-Q. 


-2Nu . c os v- 2Nv , s inv+N 


I  cosy  sinyllu 


cosv  sinv  II  v 


(M  +Q1-A‘Q2)/(2M) 
(N2+Q1-B2Q3)/(2N)_ 


(3.4) 


Since  sin(v-y)  «  sinvcosy  -  cosvsiny,  we  have 


[cosy  siny* 

-1 

_ 

»  (sin(v-u))”1 

sinv 

-smyl 

[cosv  sinv_ 

-cosv 

cosyj 

so  the  solution  of  (3.4)  is 


S' 

m 

_vl. 


smv  -Biny 


(sin(v-y)) 


'cosv  cosy 


iny  (M2+Qj-j 
2 

osyj  L(N  +QX- 


-A  Q2)/ (2m) 


BaQ3)/(2N)J 


(3.5) 


which  is  always  valid  since  sin(v-y)  *  0  because  the  circle  centres  are  not 
collinear  (this  follows  the  definitions  of  v  and  y) . 

Now  define 

a  -  (p-q-r ) /  (2/qr )  ;  6  -  (r+Q.-A2Q  )/(2»/r)  ; 

1  ^  (36) 

y  -  (q+Q1-B2Q3)/(2/i)  ; 


using  positive  square  roots  in  the  denominators  and  noting  that  p,q,r  are 
necessarily  positive.  Then  from  (3.3) 


-cos(v-y) 


-cosvcosy  -  sinvsiny 


6  -  (M2+Q1-A2Q2)/(2M) 


Y  -  (N2+Q1-B2Q3)/(2N) 


(3.7) 


Using  (j).7),  (3.5)  becomes 

£  *jj  "  (sin(v-y) ) 


jT  sinv  -siny'  “6" 

|_-cosv  cosy_  _Y. 


(3.8) 


But  from  this  and  (3.2)  we  have 


11 + vi  •  w  u;j 


(where  T  denotes  the  matrix 
transpose) 


•  Eft 


sinv  -cosv 


sinv  -sinwlfS' 


(sin(v-y)) 


sinM  cosy. 


(sin(v-y)) 


E  I] 


-cosv  cosyll  y 


(using  (3.7)) 


-  (sin(v-y))"2(S2+Y2+2a6Y) 

-  (62+Y2+2aBY)/(l-a2) 


(using  (3.7)  again) 


Substituting  for  ct,0,Y  from  (3.6)  and  simplifying  gives 

pQ2  +  qA4Q2  +  rB4Q2  ♦  (r-p-q)  A^Qj  +  (q-p-r)B2QlQ3  ♦  (p-q-r)A2B2Q2Q3 


2  2 

+  p(p-q-r)QL  +  q(q-p-r)A  Q2  +  r(r-p-q)B  Q3  +  pqr  »  0 


(3.9) 


Connecting  at  least  nine  loads  known  only  to  have  different  VRC  and 
observing  the  power  ratios  corresponding  to  each  provides  (on  dividing 
(3.9)  throughout  by  (pqr))  at  least  nine  simultaneous  inhomogeneous  linear 
equations  in  the  nine  unknown  coefficients 


Xj  *  p/(pqr)  -  1/ (qr) 

X2  -  qA4/ (pqr)  -  A4/ (pr) 
X3  -  rB4/ (pqr)  -  B4/(pq) 


X4  -  A  (r-p-q)/(pqr) 

X5  -  B2(q-p-r)/(pqr) 

X6  -  A2B2 (p-q-r ) / (pqr ) 


X?  -  (p (p-q-r)) /(pqr)  ■  (p-q-r) /(qr) 

Xg  -  (qA2 (q-p-r ) ) / (pqr)  ■  A2(q-p-r)/(pr) 
X9  ■  (rB2 (r-p-q) ) / (pqr)  -  B2 (r-p-q) / (pq) 


(3.10) 


where,  from  (3.9)  and  (3.10) 


The  set  of  equations  (3.11)  may  be  solved  by  a  standard  least-squares  method 
(see  section  6)  and  from  equation  (3.10)  can  be  deduced  the  relations 


r  ■ 


2X5  -  X?X9 

ZX^g-XjXy  '» 


r+q+(X?/X1)  ;  A 


2X4  -  X7Xg 
2X1X8-X4X7 


+*/prX^  ;  B  -  +*/pqX3 


(3.12) 


2  2 

in  which  use  has  been  made  of  A  and  B  being  positive. 


Equations  (3.11)  and  (3.12)  enable  approximations  to  the  five  required 
real  numbers  to  be  found.  The  solutions  are  likely  to  be  only  approximate 
because,  as  Engen  remarks  (6),  they  represent  conditions  of  tangency  of 
three  planes  to  the  ellipsoidal  paraboloid  represented  by  equation  (3.9) 
which  is  a  quartic  in  |w)J.  The  solutions  are  refined  by  using  a  standard 
Gauss-Newton  iterative  solution  for  the  five  unknowns  of  the  nine  (or  more) 
equations  (3.9).  (We  remark  that  since  these  five  reals  characterise  the 
ref lectometer  and  are  found  by  using  at  least  nine  DUT  known  only  to  have 
different  VRC,  they  provide  a  means  of  monitoring  the  long-term  stability 
of  the  instrument  with  loads  requiring  stability  only  during  observations 
of  the  corresponding  power  ratios  Qfc. ) 

The  purpose  of  finding  these  five  reals  p.q.r.A^.B^  is  to  enable  a 
unique  vector  w  to  be  calculated  from  the  three  observed  for  any  subse¬ 
quently  connected  load,  of  course.  This  can  be  done  by  writing 


or 


w  *  u+jv  s  w^cosy  -  jsiny)  ■  (u^+jvj) (cosy  -  jsiny) 


ar  cosu  sinuTTu^n 
|_-sinw  cospJLVjJ 

[cosy  siny*l  -P  sinv  -sinyipe*] 

(sin(v-y))'1 

-siny  cosyj  l_-cosv  cosuJLyJ 


or  from  (3.7) 


from  (3.8) 


-  [ 1  "  m 

L-cot(v-y)  csc(v-y)JLYJ 


u  •  8 


v  ■  (o$+y)/± 


/l-a2 


(3.13) 


But  from  the  definition  of  w^  we  have 


w  ■  Wj (cosy-  jsiny)  »  (dr+e)/(cr+l) 


which  is  identical  to  equation  (2.7)  derived  for  a  four-port  ref lectometer 
in  which  d  and  e  are  new  constants  of  the  six-port  (given  by 
d  »  dj (cosy  -  jsiny)  and  e  ■  ej (cosy  -  jsiny)) .  Since  a,B,Y  are  known  from 
(3.6)  for  any  Q^,  we  can  find  the  vector  v  corresponding  to  any  load  VRC, 
apart  from  the  sign  ambiguity  in  (3.13).  This  sign  ambiguity  in  the  six- 
port  ref lectometer  arises  from  that  of  the  sign  of  sin(v-y)  and  is  equivalent 
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to  not  knowing  in  advance  which  of  the  two  scalar  indications  of  the  receiver 
(or  vector  voltmeter)  of  a  four-port  ref lectooeter  represents  u  and  which 
represents  v.  A  way  of  resolving  this  by  a  slight  modification  of  the  four- 
port  calibration  method  of  section  2.2,  and  a  practical  implementation  of  a 
suitable  six-port  calibration  scheme>are  both  treated  in  the  next  section. 

3.2  SIX-PORT  CALIBRATION  USING  "3J"  STANDARDS 

With  the  aim  of  reducing  the  number  of  sources  of  uncertainty  of 
measurement  by  minimising  the  number  of  precisely-known  standards  needed, 
we  can  employ  the  algorithm  of  section  3.1  to  calibrate  a  six-port  reflecto- 
meter  with  three  such  standards  plus  a  load  whose  VRC  is  only  approximately 
known  (to  resolve  the  ambiguity  referred  to).  The  power  ratios  (k* 1,2,3) 
are  observed  when  these  four  loads  having  VRC  of  r,I”  ,r",r"'  are  connected  in 
turn  and  the  corresponding  w,w'  jw’^w"'  are  calculated  using  the  same  sign 
for  each  in  equation  (3.13).  This  determines  either  the  four  w  or  their 
complex  conjugates  w*.  Now  it  can  be  verified  by  substitution  in  equation 
(2.7)  (or  deduced  from  the  theory  of  cross-ratio  groups)  that 


(w-w”')  (w'-w"  ) 
(w-w' )  (w"-wm  ) 


(r-r”)  (r'-r" ) 

(f-F)  “(rrT-fmy 


(3.16) 


But  the  sign  of  the  imaginary  part  of  the  RHS  of  equation  (3.16)  is  known 
from  our  knowledge  of  the  four  VRC  and  if  it  disagrees  with  that  of  the  LHS 
then  the  sign  of  v  selected  in  equation  (3.13)  must  be  changed.  The  method 
will  fail  if  all  four  VRC  have  the  same  magnitude  because  the  RHS  of 
equation  (3.16)  would  then  be  real,  but  we  describe  below  a  calibration 
method  used  at  RSRE  to  avoid  this  failure.  It  also  avoids  the  apparent 
need  to  make  at  least  nine  physical  connections  to  the  ref lectometer  to 
provide  data  for  the  six-  to  four-port  reduction  algorithm. 


This  procedure  has  been  employed  to  calibrate  different  six-port 
ref lectometer s  at  wavelengths  ranging  from  metres  to  millimetres  and  it 
relies  on  a  short-  or  open-circuit  uncalibrated  variable  attenuator  to 
provide  six  or  more  different  VRC  for  use  with  the  algorithm.  Power  ratios 
Qk  are  observed  for  each  of  the  following: 


a.  Three  known  standards  connected  in  turn  to  provide  three  precisely 
known  different  VRC.  (In  millimetric  waveguide  these  have  been  an 
electroformed  short-circuit  first  connected  directly  and  then  via  two 
known  lengths  of  precision  waveguide  (spacers)  in  succession.  At  VHF 
a  short-circuit  and  two  calculable  coaxial  capacitors  have  been  used.) 


b.  A  nominally  matched  load. 


c.  A  variable  attenuator  set  to  at  least  three  settings  first  when 
unterminated  and  then  when  terminated  with  a  short-circuit.  (This 
provides  two  sets  of  three  different  VRC  for  which  the  phase  of  the 
two  sets  differs  by  more  than  90°. ) 


This  provides  a  set  of  ten  loads  known  to  have  different  VRC  and  the 
corresponding  power  ratios  are  used  in  the  algorithm  of  section  3.1  to 
calculate  A^,B^,p,q,r.  Thus  equations  (3.9) , (3.11) , (3.12)  are  used  in 
sequence  to  yield  approximate  values  and  then  successively  better  approxi¬ 
mations  are  found  by  applying  the  Gauss-Newton  technique  (outlined  in 
section  6)  to  equation  (3.9).  The  power  ratios  corresponding  to  the  first 


four  loads  (steps  a  and  b  above)  are  then  used  in  equations  (3.6)  and  (3.13) 
to  calculate  the  LHS  of  equation  (3.16),  the  known  values  of  r  being  used  in 
the  RHS  to  resolve  the  ambiguity  of  sign  of  v  in  (3.13).  The  three  precisely 
known  T  and  the  corresponding  values  of  w  are  then  used  in  equation  (2.9)  to 
complete  the  calibration  of  the  equivalent  four-port.  Thereafter,  the  VRC 
of  any  subsequently  connected  load  is  found  by  applying  equations  (3.6), 
(3.13)  and  (2.10)  in  sequence. 

In  practice,  the  nominally  matched  load  at  (b)  in  this  procedure  can 
be  provided  by  setting  the  attenuator  to  provide  its  maximum  attenuation, 
so  that  the  procedure  involves  staking  only  four  physical  connections  to  the 
measurement  port.  We  have  found  it  necessary  to  use  the  set  of  ten  VRC 
known  to  be  different  (rather  than  the  minimum  of  nine)  to  ensure  reliable 
convergence  of  the  Gauss-Newton  routine  used  in  solving  equation  (3.9).  We 
have  found  it  sufficient  to  terminate  the  process  of  successive  approximation 
of  the  five  real  numbers  A^,B^,p,q,r  when  a  precision  of  1  in  10°  has  been 
reached  for  each. 

For  a  well-designed  six-port  ref lectometer  junction,  the  squared 
distances  p,q,r  between  the  circle  centres  in  the  wj  plane  are  all  of  the 
same  order.  The  reliability  of  convergence  obtained  with  the  procedure 
described  has  been  tested  at  RSRE  by  employing  it  successfully  for  a  six- 
port  junction  maladjusted  to  the  extent  that  p  was  two  orders  less  in 
magnitude  than  q  or  r. 

As  remarked  earlier,  the  five  real  number s  found  by  this  algorithm 
provide  a  means  of  monitoring  the  stability  of  the  instrument  (since  their 
values  do  not  depend  on  known  T).  We  also  remark  that  the  redundancy 
provided  by  using  the  algorithm  minimises  the  variance  of,  in  effect,  the  w 
from  their  best  fit  values;  this  tends  to  average  out  the  effects  of  any 
departure  of  detector  response  from  linear  power  operation.  This  means 
that  the  values  of  T  which  can  be  smasured  for  loads  used  in  the  process 
(after  reduction  and  calibration  have  been  completed)  may  be  more  precise 
than  those  obtained  by  subsequent  swasuresmnts . 

3.3  CALIBRATION  IN  TERMS  OF  A  COMMON  Z 

o 

We  have  so  far  made  the  assumption  in  deriving  calibration  procedures 
that  there  is  no  difficulty  in  realising  a  set  of  practical  calibration 
standards  to  provide  at  least  three  precisely-known  different  VRC  norma¬ 
lised  to  the  same  ZQ.  For  example,  in  calibrating  six-ports  with  a  wave¬ 
guide  short-circuit  and  two  known  lengths  of  guide  we  have  approximated  to 
truth  by  considering  the  cross-sections  of  the  waveguide  to  be  identical 
and  the  standards  to  be  lossless.  The  precision  of  measurement  achieved  at 
RSRE  with  six-port  ref lactometers  is  greater  than  that  of  these  approxima¬ 
tions  for  coomercia 1-quality  waveguide  and  the  performance  recently  reported 
of  an  improved  Sttvable  short-circuit  (18)  suggests  that  it  is  worthwhile 
resKiving  at  least  the  need  to  assume  lossless  waveguide  (although  the  VRC 
of  the  short-circuit  must  be  asstaned  known).  This  can  be  done  by  applying 
to  a  six-port  a  method  developed  by  da  Silva  and  McPhun  of  calibrating  a 
four-port  ref lectometer  (14). 

The  method  assumes  that  the  VRC  of  a  reflector  (the  short-circuit) 
retains  a  constant  value  rQ  irrespective  of  its  position  along  a  uniform 
transmission  line  or  waveguide.  If  the  power  ratios  Qfc  are  observed  at 
scale  settings  of  * •  0,  (£0+i),  (t0+2£),  (£0+3£)  along  the  line  measured 
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from  Che  connection,  Chen  Che  VRC  at  Che  connection  will  be  rp, 
r'  ■  r0  exp(-2 61),  T"  •  rQ  exp(-461),  r'"  ■  r0exp(-66i),  6  being  the  propa¬ 
gation  constant  of  the  line.  Corresponding  vector  indications  wotv' 
can  be  calculated  from  the  power  ratios  and,  from  equation  (3.16) 


(wo-w"')  (w'-w"  ) 

(w°-w*)(w»-w,n) 

o 


(r^Dfr'-r" ) 
TT^TTF^PT 

O 


(3.17) 


Writing  x  ■  exp (-261)  then  the  RHS  of  this  equation  is  (x^+x+D/x  and 
denoting  the  LHS  (calculated  using  the  positive  root  in  (3.13))  by  W,  then 


•  (x*+x+l)/x  and  so 


(W-l±  /(W-l)  -4)/2 


(3.18) 


The  two  roots  of  (3.18)  are  reciprocals  of  each  other  and  only  the  one  with 
modulus  less  than  unity  is  admissible,  since  the  transmission  line  attenu¬ 
ates.  The  other  possible  value  of  the  LHS  of  (3.17)  is  the  complex  conju¬ 
gate  W*  and  the  version  of  (3.18)  corresponding  to  this  has  the  admissible 
root  x* ,  where  x  is  the  admissible  root  obtained  fromW.  To  choose  between 
x  and  x*  we  remember  that  theory  predicts  that  for  smooth  conductors 
^  a  (jk),  where  k  is  the  free-medium  angular  wave  number  in  the  transmission 
line  (usually  written,  with  symbol  usage  different  to  that  of  this  paper, 
bt/v  or  w/eu/c) .  Let  Log  represent  the  complex  logarithm  function  with 
phase  0  mapped  into  the  range  -»  <  8  <  ir,  then  if  x,  rather  than  x*,  is 
correct  we  must  have 


Log  x  *  -2jki  +  2m 
where  n  is  a  suitable  integer. 


(3.19) 


In  particular  the  signs  of  the  imaginary  parts  of  the  LHS  and  RHS  of 
(3.17)  will  agree.  If  they  do,  we  were  correct  in  choosing  x,  W  and  the 
positive  sign  in  (3.13);  if  not  then  we  should  have  chosen  x*,  V*  and  the 
negative  sign  in  (3.13).  This  is  sufficient  to  resolve  the  ambiguity,  for 
the  sign  of  (3.13)  is  independent  of  T;  henceforth  we  shall  denote  the 
correct  root  by  x. 

We  now  know  x,  so 

2  3 

r  ,  r '  ■  r  x,  r"  ■  r  x  r'"  ■  r  x 

o  o  o'  o 

are  all  known  multiples  of  r0,  and  we  know  the  corresponding  w0,w,,w",w"’ 
so  that,  from  equation  (2.7) 

dr  ♦  e  -  w  cr  ■  w  ;  dTx  +  e-  w'cr  x  ■  w'  ; 
o  ooo  o  o 

dr  x2  ♦  e  -w”cr  x^  -  w"  ;  dr  x^  ♦  e  -  w"'cr  x^  -  w'" 
o  o  o  o 


which  are  four  simultaneous  linear  equations  in  the  three  unknowns  (cro), 
(dr0),  e  which  can  be  solved  by  a  standard  least-squares  method.  This 
knowledge  allows  the  ratio  of  any  two  VRC  r,  r’  to  be  calculated  from  the 
corresponding  indications  w,  w'  obtained  from  the  observed  power  ratios 
for,  from  equation  (2.7), 


TCTTr  at  'i ! 


r  a  W~S 
r  (dr  )-w(cr  ) 


But 


(r/r  )/(r’/r  ) 

o  o 


so  that  the  ratio  can  be  found, 
must  be  known. 


To  obtain  any  T  absolutely,  however,  r£ 


Variations  of  results  obtained  in  experiments  at  RSRE  using  a  sliding 
short-circuit  to  implement  this  calibration  method  using  7  mm  bore  precision 
50  ohm  airline  as  a  standard  showed  that  the  assumption  of  a  constant 
ro  *  -1  for  all  short-circuit  positions  was  inadequate.  Nevertheless, 
apart  from  its  possible  application  were  improved  short-circuits  available, 
it  is  of  interest  because  of  its  affinity  with  the  method  of  calibrating  a 
dual  six-port  network  analyser  described  later.  This  uses  as  the  standard 
a  known  length  £  of  uniform  line  (or  waveguide);  a  short-circuit  which 
is  assumed  only  to  have  the  same  VRC  for  each  connection  is  also  used  but 
its  r?  is  measured  by  the  calibration  process,  as  is  the  x  *  exp(-22£)  of 
the  line. 

3.4  POWER  FLUX  FROM  THE  MEASUREMENT  PORT 

So  far  we  have  assumed  that  each  of  the  four  detectors  provides  an 
indication  only  proportional  to  the  power  absorbed  by  it.  If  now  we  have 
a  power  meter  (ie  one  which  correctly  indicates  the  power  absorbed)  to 
serve  as  a  standard  of  power  then  a  six-port  reflectometer  already  calibrated 
for  VRC  measurement  can  be  calibrated  to  indicate  the  power  flux  emerging 
from  its  measurement  port.  It  may  then  be  used  for  comparing  power  meters 
by  connecting  them  in  turn,  eg  for  comparing  the  effective  efficiencies  of 
boloaieter  mounts. 


«3r  .  e3|2|b2|2 


!B3lVcr|2|i>2|2 


If  the  power  absorbed  by  the  standard  power  meter  having  VRC  rL  is  PL  then 
PL  ' 

■  =o(1  ’  lrJ‘)lb2l' 


'2  -  l*2|2) 


(3.20) 


where  GQ  is  the  conductive  component  of  the  admittance  characterising  the 
mode  transmitted  at  the  measurement  port. 


By  division 
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mtu  *7  *s  7*  ^  ■?  ■?  i\  ■-■.  :■.  '■r.v 


, 


4 


P 


L 


°o  »-|rL|2) 
|8j|2  |i*crL|2 


(3.21) 


Substituting  for  from  equation  (2.7) 


Writing 


K 


G 

o 


|6jI2U-c.|2 


as  a  constant  characterising  the  ref lectometer  gives 

PL  -  K(|d-wc|2  -  |v-e|2)P3  (3.22) 

Both  K  in  (3.22)  or  G0/ | 83 | 2  in  (3.21)  characterise  the  ref lectometer  and 
can  be  found  by  observing  and  the  indication  of  the  standard  power  meter 
when  connected  to  the  measurement  port;  thereafter  the  power  absorbed  by  any 
DUT  can  be  found  from  either  of  equations  (3.21)  or  (3.22). 


THE  DUAL  REFLECTOMETER  NETWORK  ANALYSER 


4.1  PRINCIPLES,  CALIBRATION  PROCEDURE  AND  NOTATION 


So  far  we  have  been  concerned  only  with  measuring  the  VRC  and  power 
absorbed  by  one-port  DUTs,  but  the  scattering  coefficients  of  two-port 
junctions  can  be  measured  with  two  ref lectometer s  connected  to  the  same 
source  and  a  phase-changer,  as  shown  in  Figure  3.  We  illustrate 

the  principle  of  this  by  assuming  that  the  ref lectometer s  have  been  cali¬ 
brated  so  that  their  indications  are  apparent  VRCs  in  which  the  "reflected" 
waves  comprise  both  those  reflected  by  and  transmitted  through  the  two- 
port  junction  under  test.  Thus,  using  the  notation  of  Figure  3,  we  have 
rA  »  aA/bA  and  Tg  -  ag/bg.  Denoting  the  scattering  coefficients  of  the  DUT 
by  s£j  in  the  usual  way  gives 


*A  *  *llbA  *  *12bB  1 
*B  "  *21bA  +  *22bB  ) 


*°  tb,t  rA  '  *11  *  *12  K 


,ni  rB  ’  *21  +  *22 


Eliminating  (bA/bg)  from  these  two  equations  gives 


rB*ll  *  rA*22  * 


(,12*2l",ll*22) 


r  r 
Vb 


(4.2) 
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By  setting  the  adjustable  phase-shifter  to  three  (or  more)  values  known 
only  to  provide  different  phase  relations  between  and  bg,  there  are 
obtained  three  (or  more)  sets  of  (rg)i,  (rg)j,  "  1*2.3  ..)  which  result 
in  a  set  of  three  (or  more)  simultaneous  equations  (4.2)  which  can  be 
solved  for  sn,  822  and  the  product  812*21*  This  allows  determination  of 
the  scattering  coefficients  of  any  two-port  junction  that  is  reciprocal  (in 
the  sense  that  *12  “  *21  w*ien  both  ports  transmit  the  same  mode  (15)), 
provided  we  ignore  for  the  moment  the  ambiguity  of  phase  of  *12*  We  remark 
that  the  phase-shifter  needs  constancy  only  during  each  measurement  and 
does  not  need  to  accurately  repeat  its  settings  for  a  sequence  of  measure¬ 
ments  (unless,  as  we  shall  see,  there  is  a  requirement  to  measure  phase 
angles  £812  *nd  £*21  *  non-reciprocal  two-port  junction).  We  also  remark 

that  the  reflectometers  need  not  be  six-ports  and  that  a  phase-shifter  with 
four  settings  allows  the  measured  result  to  be  the  best  fit  to  four  sets  of 
power  ratio  data  (since  it  is  a  least-squares  solution  of  an  overdetermined 
set  of  equations  (4.2)). 

Each  apparent  VRC  indicated  during  measurement  of  a  passive,  two-port 
junction  stay  be  greater  than  unity,  as  can  be  seen  by  considering  the  "null 
two-port  junction"  (ie  when  the  two  reflectometers  are  connected  together) 
for  then  a^  and  bg  are  equal  (r^  ■  l/r-l.  This  implies  that  the  reflecto- 
meter  calibration  methods  so  far  derived  may  be  inappropriate  because  they 
employ  standards  of  known  VRC  of  j  r {  <1.  Recalling  the  coranents  following 
equation  (2.10),  this  means  that  either  the  excitation  of  each  reflectometer 
must  be  halved  after  calibration  to  ensure  that  the  detectors  always  operate 
in  their  linear  range  or  it  must  be  assumed  that  their  response  remains 
linear  when  |f|  >  1.  The  calibration  procedure  we  now  consider  is  essentially 
the  one  described  in  (8)  and  utilises  connection  of  both  reflectometer 
measurement  ports  to  each  other  and  calibrates  both  instruments  in  terms  of 
just  one  precisely  known  standard,  namely  a  known  length  i  of  precision 
waveguide  or  air-dielectric  transmission  line  having  the  same  uniform  cross- 


section  as  the  measurement  ports.  The  steps  in  the  calibration  procedure  are 
illustrated  in  Figure  4  and  they  consist  of  observing  the  power  ratios: 

a.  when  the  two  measurement  ports  are  connected  together  (condition  1) 

b.  when  the  two  measurement  ports  are  connected  via  the  length  l  of 
standard  waveguide  or  coaxial  line  (condition  2); 

c.  when  they  are  connected  via  a  stable  but  unknown  two-port,  such 
as  a  6  dB  or  10  dB  pad; 

d.  when  there  is  zero  transmission  between  the  measurement  ports 
(condition  3);  this  is  achieved  by  terminating  both  ports  with  different 
unknown  reflecting  terminations,  which  are  exchanged  to  allow  two  sets 

of  observations  to  be  made. 


Refl’r 

Refl’r 

A 

B 

(a)  measurement  ports 
-  connected  directly 
together  ( condition  1 ) 


Refl’r 

Standard 

Refl*  r 

A 

'  line  5 

B 

(b)  connected  via  standard 
line  (condition  2) 


Refl'r  _Jtw31  ,  Refl’r 
A  “Hport  B 


(c)  connected  via  a 
reciprocal  two -port 


Refl’r 

Load  1 

.  Refl’r 

A 

Load  2 

B 

(d)  both  reflectometers 
terminated:  measured 
also  with  loads  1  &2 
interchanged  (condition  3) 


FIG.  4.  CALIBRATION  PROCEDURE  FOR  DUAL  -  REFLECTOMETER 


NETWORK  ANALYSER 


Only  the  three  conditions  numbered  1  to  3  in  brackets  are  essential  to  the 
calibration  process  (except  for  measurement  of  the  phase  angles  £Si2  and  £S2i 
of  non-reciprocal  two-port  junctions).  The  object  of  using  additional  unknown 
one-  and  two-port  devices  is  to  provide  overdetermined  sets  of  equations  to 
enhance  the  measurement  precision.  Every  least-squares  solution  implies  a 
consistency  check  and  if  the  unknowns  are  stable  they  also  provide  checks  of 
consistency  of  calibration. 


All  four  states  of  phase  excitation  are  used  during  steps  a,  b,  c , 
but  only  one  is  needed  when  there  is  zero  transmission  at  step  d.  There 
are  therefore  fourteen  sets  of  observed  power  ratios  for  each  six-port 
available  as  data  for  applying  the  algorithm  of  section  3.1  to  reduce  each 
to  an  equivalent  four-port  ref lectometer  (complete  with  its  own  ambiguity 
of  sign  of  v  resulting  from  equation  (3.13)).  These  equivalent  four-ports 
are  modelled  by  the  two-port  error  box  representation  of  section  2.2  and 
the  notation  illustrated  in  Figure  5  is  used  with  this  representation. 

The  waves  incident  on  and  emergent  from  the  notional  perfect  ref lectometer s 
A  and  B  are  denoted  by  a^,  b^  and  as,  bg;  the  vector  indications  of  these 
notional  ref lectometers  by  w^  *  aA^A  an<*  WB  “  aB^B»  an<*  rA»  rB  are  used 
for  the  VRCs  at  the  true  measurement  ports.  With  this  notation,  (4.1)  can 
be  applied  to  a  two-port  junction  connected  between  the  two  notional 
measurement  planes  to  relate  wA,  wg  to  its  scattering  coefficients  (by 
analogy  with  the  derivation  of  (4.2),  giving 

Vll  *  V22  *  U12'21 '  sll*22)  '  WAVB  (4'3) 


During  measurement  a  DUT  can  be  considered  as  being  in  cascade  between  the 
two  error  boxes  so  it  is  convenient  to  express  (4.1)  in  a  different  matrix 
form 


M 

1 

(81282rSll822) 

*11 

Y 

1 

JL 

S21 

-s22 

1 

m 

8B. 

(4.4) 


The  square  matrix  and  its  scalar  multiplier  in  equation  (4.4)  is  the  so- 
called  cascading  matrix  (15)  and  this  has  the  important  property  that  the 
cascading  matrix  of  two  (or  more)  two-port  junctions  connected  in  cascade 
is  the  ordinary  matrix  product  of  their  individual  cascading  matrices. 


We  now  need  to  introduce  further  notation  (illustrated  in  Figure  5)  to 
handle  the  observations  made  during  conditions  1  to  3  of  the  calibration 
procedure. 


“Perfect"  Error 
ref  lectometer  box 


Error  “Perfect* 
box  reflectometer 


aR  and  R 
A  B 


are  cascading  matrices  of  error  boxes 


FIGS.  NOTATION  USED  FOR  DUAL  SIX -PORT  NETWORK  ANALYSER 
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Denoting  the  cascading  matrices  for  error  boxes  A  and  B  by  (ar),(br)  * 
respectively,  and  that  of  the  standard  line  by  (LR),  and  using  superscripts 
1,2,3  to  represent  for  these  three  conditions  the  vector  indications  v^.vg, 
enables  us  to  write  Iwg  for  these  indications  when  the  two  reflecto- 

meters  are  connected  together  (ie  connected  via  the  "null  two-port"),  and 
V  ^wg  when  they  are  connected  via  the  standard  length  of  line,  etc. 

This  notation,  plus  the  observation  that  the  cascading  matrix  of  the  null 
two-port  is  just  the  (2x2)  identity  matrix  and  can  therefore  be  suppressed. 


will  enable  us  to  set  up  cascading  matrices  to  relate  ^w^, 

(AR>,  W.  (lr)« 


^B*  2wA»  ^  to 


The  two  six-port  ref lectometers  have  eight  power  meters  between  them 
and  these  give  seven  independent  power  ratios,  but  the  two  sets  of  Qg  (as 
defined  in  section  3.1)  account  for  only  six  of  them.  For  measurements  of 
non-reciprocal  two-port  junctions  we  also  need  the  seventh  and  the  defini¬ 
tions  of  the  power  ratios  Qg  must  be  extended  to  include  it.  If  the  indi¬ 
cation  from  detector  D3  of  Figure  1  is  hereafter  called  Pg  (instead  of  P3) 
and  the  others  renumbered  1,2,3  for  each  six-port  then  we  have,  in  an 
obvious  notation. 


QAk 

P  /P 

Ak7  AR 

(k  - 

1,2,3)  ’ 

| 

QBk 

P  /P 

Bk'  BR 

(k  - 

1,2,3) 

►  (4.5) 

QAB  * 

P  /P 
BR/rAR 

Note  that  a  measurement  set  of  ratios  includes  only  one  Q^g. 

Having  illustrated  the  principle  of  operation,  described  our  calibra¬ 
tion  procedure,  and  established  our  notation,  we  can  proceed  to  derive  our 
calibration  theory.  We  tackle  first  the  not  inconsiderable  problem  of 
resolving  the  ambiguities  of  sign  of  v  for  the  two  equivalent  four-port 
ref lectometers. 


4.2  RESOLVING  SIGN  AMBIGUITIES  FOR  DUAL  SIX-PORT  NETWORK  ANALYSER 


Having  reduced  each  six-port  to  an  equivalent  four-port  ref lectometer, 
we  apply  equation  (3.16)  to  the  measurements  made  with  condition  1  of  the 
calibration  procedure  to  give 


,1  1  W1  1  . 

(  'I*'  V  (  ”2A*  W3A 

TI  I  77l  I  7 

<  ”1A'  W2A  (  W3A~  V 


(lriA-lr<,A)(lr2A-lr3A) 

<lriA-lr2A)(lr3A-lr4A> 


[  riB~  r4B^(  r2B~  1W 
<lrilflr2B)(lr3lflr4B> 


,1  l  ,,1  1  , 

(  “IB"  U4B)(  W2B~  V 

71  i  7a  I  7 

<  “IB-  W2B  (  W3B-  ”4B> 


\ 


(4.6) 


because  for  the  null  two-port  we  have  *riB  *  1  (i  ■  1,2, 3, 4)  for  all 

four  phase  states.  The  signs  of  the  imaginary  parts  of  the  cross-ratios  in 
W'x  and  Wg  are  compared  and  if  they  disagree  then  the  negative  sign  is  taken 
in  equation  (3.13)  for  either  w^  or  wg  (at  this  stage  it  does  not  matter 
which).  This  does  not  guarantee  that  both  signs  are  now  right  but  it  does 
guarantee  that  either  they  are  both  right  or  both  are  wrong;  that  is, we  are 

now  using  consistently  for  A  and  B  either  w  or  w*  (but  we  do  not  yet  know 

which).  We  assume  for  the  moment  that  we  have  w  and  correct  later  if  this 
is  found  to  be  wrong.  It  is  worth  noting  that  (4.6),  unlike  the  use  of 

(3.16)  in  section  3.2,  is  exact  so  that  (after  any  necessary  sign  change) 

the  two  cross-ratios  in  w  in  (4.6)  should  be  exactly  equal.  This  is  one 
of  many  consistency  checks  that  this  calibration  method  automatically 
provides. 

Next,  we  calculate  the  elements  of  the  cascading  matrices  of  the  two- 
ports  between  the  notional  measurement  planes  of  w^  and  wg,  first  when  the 
true  measurement  planes  are  connected  together  (condition  1)  and  then  when 
they  are  connected  via  the  length  of  standard  transmission  line  (condition 
2).  To  do  this,  we  apply  (4.3)  for  each  of  the  four  measurements  (with 
different  phase  excitation)  to  each  of  the  two  conditions;  for  example  for 
condition  1,  this  gives 


11  A  1  1  A  ,1  1  11.11 

■a  *11  *  ”iA  "22  *  <  *12  s21*  'll  *22)  '  wiA  wiB 


(4.7) 


(i  -  1,2, 3, 4) 

which  provides  four  simultaneous  inhomogeneous  linear  equations  in  the  thre« 
unknowns  Is^,  *822  and  (l8i2ls21"  ^8llls22^;  these  unknowns  may  now  be 
found  by  a  standard  method.  From  (4.4),  we  immediately  have  the  elements 
of  the  corresponding  cascading  matrix,  except  for  a  common  factor  (l/1S2i); 
this  is  ignored  because  (as  we  shall  see)  only  the  ratios  of  the  elements 
are  significant  provided  ^821  *  0  (which  is  obviously  true  whenever  trans¬ 
mission  between  the  two  notional  planes  of  w^  and  wg  is  possible).  Applying 
the  same  procedure  to  the  measurements  associated  with  condition  2  gives  the 
elements  of  its  corresponding  cascading  matrix,  apart  from  an  ignorable 
factor  of  (l/^82i).  (We  remark  that,  since  the  determinant  of  the  elements 
of  the  cascading  matrix  of  a  reciprocal  two-port  junction  with  both  ports 
referred  to  the  same  reference  ZQ  is  unity,  the  elements  can  be  determined 
absolutely,  apart  from  a  sign  ambiguity;  however,  the  assumption  of 
reciprocity  is  an  unnecessary  restriction  and  yet  another  sign  ambiguity  is 
obviously  undesirable.) 

Denoting  the  cascading  matrices  just  found  for  calibration  conditions  1 
and  2  by  *R  and  respectively,  we  see  that  they  satisfy  the  equations 

lR  "  ls21(AR)(BR)  and  2r  “  2®21(AR)(LR)(BR)  (4,8) 


where  the  cascading  matrix  of  the  null  two-port  has  been  suppressed.  We 
also  have,  using  the  notation  of  section  3.3  in  (4.4) 


Now  define  a  matrix  T  s  2R(1R)”1  det(1R).  (As  remarked  earlier,  det(*R) 
is  usually  unity  but  we  do  not  need  to  assume  that  it  is,  only  that  it  is 
*  0.  For  this,  the  sufficient  condition  is  the  physically  essential  relation 
^s21  ^  0‘)  The  elements  of  T  can  be  calculated  from  the  known  elements  of 
1r  and  2R;  using  ij  to  denote  elements,  we  obtain 


2  1  2  1 
Tij  *  Rij  R(3-j) (3-j)  "  Ri(3-j)  R(3-j) j 


(4.10) 


But  we  also  have,  on  eliminating  gR  between  the  two  equations  (4.8), 


2s21<AR)(LR)((1/1s21)(AR)"1  1r) 


so 


T  -  det(1R)2R(1R)_1  -  det(XR) (2s21/ls21) (AR) (jR)(aR) 

or,  on  writing 

det(1R) (2s2^/^s2i) 


as  <p  (a  number,  not  a  matrix), 
T(AR)  “  *(AR)(LR) 


Finally,  we  multiply  this  out,  using  the  suffix  notation  and 
ting  from  (4.9)  for  LRij »  to  give 

substitu- 

T11  AR11  +  T12  AR21 

- 

*ARll 

exp (-62.) 

(4.11.1) 

T21  AR11  +  T22  AR21 

m 

*AR21 

exp(-6£) 

(4.11.2) 

T11  AR12  +  T12  AR22 

- 

*AR12 

exp (6£) 

(4.11.3) 

T2l  AR12  +  T22  AR22 

- 

^AR22 

exp(6£) 

(4.11.4) 

Dividing  (4.11.1)  by  ARn  and  (4.11.2)  by  ar21  giye8 

T11  +  T12(AR21/AR11)  “  *  exP("62'>  “  T21<AR11/AR21)  +  T22 


T2i(ARll/AR21)  +  (T22"T11) (AR11/AR21) 


(4.12.1) 


Treating  (4.11.3)  and  (4.11.4)  similarly  gives 


T21 ^AR12^AR22^  +  ^T22~^ll^  ^AR12^AR22^ 


(4.12.2) 


Also,  on  dividing  (4.11.1)  by  (4.11.4)  and  recalling  the  notation  of 
section  3.3,  we  have 


x  -  exp (-261) 


T11*T12*AR21^AR11* 

T22+T21(AR12/AR22) 


(4.13) 


Now  (4.12.1)  and  (4.12.2)  show  that  (ARll)/(AR2l)  an<*  (ar12)/(ar22) 
satisfy  the  same  quadratic  equation.  There  are  two  possibilities  here; 
either  the  two  element  ratios  are  given  by  the  same  root  of  the  quadratic, 
or  one  of  them  is  equal  to  one  root  and  the  other  to  the  other  root.  We 
can  reject  the  first  possibility  because  it  requires  the  two  element  ratios 
to  be  equal  which,  from  (4.4)  requires  that  ^s^  As21  “  this  is  plainly 
non-physical.  Taking  the  second  possibility  then,  and  solving  (4.12)  by 
formula,  gives 


(AR12/AR22) »  ^AR1 1  ^ AR2 1  ^ 


Tll“T22±,^Tll“T22)  +4T12T21 


2T 


(4.14) 
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Substituting  these  values  into  (4.13)  gives 


exp(-26i) 


Tn,T22^(Tn-T22)2<4T 

Tn*T22l/<TirT22)2’4T 


T 

12x21 

T 

12a21 


where  the  upper  signs  in  the  numerator  and  denominator  correspond  to  taking 
the  upper  sign  in  (4.14)  for  (ar12^/^Ar22^ * 

The  two  possible  values  are  reciprocal  so,  as  in  section  3.3, 

we  reject  the  one  with  modulus  >  1.  (This  requires  a  precision  of  measure¬ 
ment  that  enables  this  discrimination  to  be  effective  for  the  few  centi¬ 
metres  of  air  dielectric  line  used  as  the  calibration  standard.  It  has  been 
found  to  be  effective  for  75  mm  of  7  mm  bore  precision  coaxial  line  giving 
about  0.004  dB  attenuation;  alternative  methods  similar  to  that  of  section 
3.2  are  available,  however,  for  systems  with  insufficient  precision.)  Now, 
consider  the  effect  of  choosing  the  wrong  sign  in  (4.6),  so  that  w*  has  been 

used  throughout  instead  of  w.  Reviewing  the  derivation  so  far,  it  follows 

that  the  effect  would  be  to  give  *s*  and  2s*  for  *s  and  2s,  and  similarly 

*R*,  2R*,  T*,  (^Ri2/ar22)*»  ^Ar11 ^Ar21^*  an<*  x*  instead  of  their  unconju¬ 

gated  forms.  We  resolve  this  by  using  the  method  of  section  3.3  to  choose 
between  x  and  x*  and  all  other  quantities  are  automatically  resolved  as 
well. 


It  is  important  to  note,  however,  that  with  6  approximately  equal  to  jk, 
the  discrimination  depends  on  the  ability  to  distinguish  between  exp(-2jk£) 
end  (exp(-2jkl))*  -  exp(2jk£).  If  these  quantities  are  essentially  equal, 
as  will  happen  if  the  length  l  is  too  close  to  an  integral  number  of  quarter- 
wavelengths,  then  the  method  will  not  work.  The  determined  difference 
between  6  and  jk  affords  a  very  sensitive  test  of  the  repeatability  of 
calibration. 

All  the  ambiguous  signs  having  been  sorted  out  for  the  two  equivalent 
four-port  ref lactometers,  we  can  now  derive  the  theory  for  their  calibration. 
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4.3  DUAL  REFLECTOMETER  CALIBRATION  THEORY 

We  begin  by  applying  equation  (4.3)  to  transform  from  wA  to  I*A  instead 
of  from  wA  to  w-.  Prefixing  S£ j  by  A  to  represent  the  elements  of  the 
scattering  matrix  of  error  box  A  we  have,  remembering  that  we  must  invert 
rA  to  allow  for  the  reversed  measurement  direction  (as  shown  by  the  arrows 
in  Figure  3). 

Tkk*ll  +  WAAS22+  (AS12AS2l"ASll  AS22)  "  Va 

Using  (4.4),  this  becomes 

rAlAR12  "  WAAR21  +  AR11  ’  Va  AR22 


or 


WA_(AR12/AR22) 


(4.15) 


(4.16) 


A  (AR11/AR22) (1_WA(AR21/AR11)) 

Comparing  this  with  equation  (2.7),  we  obtain 

(cA/dA)  "  1/(AR11/AR21)  ;  dA  “  (AR11/AR22)  ; 

eA  "  <AR12/AR22)  ; 

From  (4.14)  we  have  now  determined  (cA/dA)  and  eA. 

Apply  (4.3)  to  transform  from  Tg  to  wB  instead  of  from  wA  to  w^;  this 
time  we  must  invert  Tg  to  allow  for  the  reversed  measurement  direction  and, 
by  the  procedure  just  used,  we  obtain 

WBBSll+  rBlBS22  +  (B*12  BS21  "BS11  BS22)  "  FB  WB 


or 


WB  BR12  ”  rB  BR21  +  BR11  *  PB  WB  BR22 


or 


B 


wb'*,(BR21/BR22) 
BR22^ 


(4.17) 


B  B  12  B  11 

which  gives,  by  comparison  with  (2.7)  (or  with  (4.15)  and  C4.16)) 
“  "(BR12/BR11)  J  dB  "  <BR11/BR22)  5 

eB  *  ■(BR21/BR22) 


(4.18) 
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To  obtain  information  about  the  quantities  in  (4.18)  we  appe 
(4.8);  substituting  from  (4.16)  and  (4.18)  into  *R  *  1*21^AR^BR^ 


appeal  to 
(bR)  gives 


*11  *4.2  .  dA  ®A 

“  *21 

lR  Hi  cl 

.  r21  k22.  .A  . 


If  we  premultiply  this  by  C^R)-1  in  the  form 


dB  ~CB 


“eB  1 


r  1  “eA- 

lR)  AR22 

L'ca  VI 


(noting  that  the  physical  necessity  that  ^$21  and  4*12  are  non-zero  means 
that  the  determinant  is  also  non-zero) ,  we  obtain 


*  *c,l  i 

I  "  (AR22(  ®21 

K  1  J 


BR22  det(AR)) 


,  r 1  -vt1ru  1ri2i 

*l) 

,  1  1 

L'CA  dAJL  *21  R22-» 


(4.19) 


where  all  the  divisions  in  the  numerical  factor  are  legitimate,  including 
that  by  gR22  which  is  (b*21)~*  similarly,  since 


AR22  “  Vn*  *  0 
the  factor  is  not  zero  either.  Writing 

6  *  AR22(  *21  BR22  det(AR)) 
into  (4.19)  then  gives 

-CB  -  8  (  *12-«A  *22>  i  dB  ■  8  (  Rn~ex  *21)  1 

-*B  *  6<-cA1ru"!A1r21)  •  1  ■  8(-8a1r12*<1A1r22) 


Dividing  out  0, 


<w 


i.  l*  lv  l» 

*12  eA  *22  .  d  _  *11  eA  *21 

lRH“eA  R21  dA(  R22“(cA/dA)  *12* 


^21*  RH 

lR22“CcA/dA)  R12 


(4.20) 
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Since  we  already  know  (c^/dA)  and  e* ,  equations  (4.20),  together  with  the 
IRi:  defined  at  (4.8),  give  us  (cg/dg),  eg  and  (dgdA)  (on  multiplying  the 
second  of  (4.20)  by  dA). 

From  (2.7),  (4.15),  (4.16),  (4.17)  and  (4.18), 

r*  ‘  'MW1  “*  r»  '  Vl-VVV>  <4-n) 

with  all  the  parameters  known  except  dA  and  dg  whose  product  is  known. 

Multiplying  equations  (4.21)  together  gives  an  expression  for  (rArg) 
which  is  completely  determined  (since  the  factors  dA  and  dg  combine  into 
their  known  product  and  everything  else  is  already  known) .  We  now  invoke 
condition  3  of  the  calibration  procedure  in  which  the  same  fixed  load  is 
presented  first  to  one  measurement  port  and  then  to  the  other.  The  FA  of 
one  of  these  observations  is  clearly  equal  to  the  Tg  of  the  other  (and  vice 
versa,  since  the  realisation  in  calibration  step  d  consists  of  two  fixed 
loads  which  are  exchanged)  and  their  product  is  Tx^,  where  rx  is  the  VRC 
of  the  relevant  termination.  A  trivial  amount  of  information  about  the 
nature  of  the  fixed  termination(s)  is  sufficient  to  determine  the  correct 
square  root  of  rx2  (eg  if  rx2  a  1  4-jo  one  can  distinguish  between  short- 
and  open-circuits  by  inspection).  Then,  using  equations  (4.21)  separately 
dA  and  dg  can  be  found  (and  measurements  on  both  fixed  terminations  provide 
for  a  check  by  comparison  of  results).  The  cA,  dA,  eA  and  eg,  dg,  eg 
characterise  the  two  ref lec tome ter s  so  that  calibration  is  now  complete. 

The  rA  or  Tg  at  the  measurement  port  of  each  six-port  ref lectometer 
can  be  found  from  the  observed  by  determining  a,  8  and  y  from  equation 
(3.6),  finding  the  appropriate  w  from  (3.13)  (since  the  appropriate  sign 
for  v  is  known),  and  finally  obtaining  T  from  (4.21).  The  TA£  and  Tgj 
obtained  for  the  settings  of  the  adjustable  phase  —  changer  can  then  be  used 
to  set  up  a  set  of  at  least  three  simultaneous  inhomogeneous  equations 
(4.2)  which  are  solved  by  standard  methods  for  su,  S22  end  the  product 
*12*21  °f  *  two-port  junction.  If  this  is  reciprocal  In  the  sense  that 
*12  “  *21  then  i*  can  be  completely  characterised  provided  plausible 
arguments,  like  those  of  section  4.2,  may  be  relied  on  to  resolve  the 
square-root  ambiguity  (but  see  sections  4.4  and  4.5). 


The  calibration  procedure  of  section  4.1  and  Figure  4  provides  for 
these  measurements  on  all  the  one-  and  two-port  junctions  used  and  the 
results  serve  to  monitor  the  consistency  of  calibration. 

4.4  CALIBRATION  THEORY  FOR  MEASURING  | si2 I  AND  | S21 I  OF  A  NON-RECIPROCAL 

TWO-PORT  JUNCTION  TRANSMITTING  THE  SAME  MODE  AT  BOTH  PORTS 

So  far,  we  have  not  made  use  of  the  seventh  power  ratio  Q^g  »  Pbr/par 
defined  at  (4.5);  this  provides  another  calibration  constant  which  enables 
| S|_2 1  and  | S21 1  of  a  two-port  junction  to  be  determined.  (Thus  the  assump¬ 
tion  of  reciprocity  8^2*21  “  1  can  be  verified  within  an  ambiguity  of  an 
even  multiple  of  s,  which  should  be  resolvable  from  physical  considerations.) 

From  equation  (3.20)  the  power  flux  from  the  measurement  port  of 
reflectometer  A  is 


pa  *  Go<HrAni»Al2 


(4.22 


But  from  (3*22)  and  the  definitions  at  (4.5) 


PA  '  VaR^a'aI*  -  I»A-a|2> 
Substituting  for  Ig  from  (4.21)  this  becomes 

pa  ■  VarIVVa^Hi'aI2) 

l*>A*2  •  VarIWa'X  <flOT  <‘-22» 


Evaluating  Pg»  the  power  flux  from  measurement  port  B,  similarly  gives 


bA 

2 

ka 

Wa 

bB 

“  *B 

Wb 

(4.23) 


But  (K&/Kg)  is  a  characteristic  of  the  network  analyser  and  can  be  found 
from  the  observations  of  the  power  ratio  Qa^  when  the  two  ref lectometers 
are  connected  together  (calibration  condition  J.)  for  then,  using  primes 
in  equation  (4.22)  to  indicate  this  condition 


Butp| 


Pg  so  that,  from  this  equation,  (4.21)  and  (4.23) 


(4.24) 


Thus  (Kg/Kp)  can  be  found  by  observing  QAg  ■  PgR/PAR  (as  in  (4.5)) 
during  calibration  condition  1.  Since  it  is  a  constant  for  the  network 
analyser,  and  since  observations  are  made  with  at  least  three  settings  of 
the  phase-changer,  it  provides  another  check  of  calibration  consistency. 
(The  usual  least-squares  solution  for  the  best  value  conveniently  reduces 
to  the  arithmetic  mean  of  the  separately  calculated  values.) 


Prom  equations  (4.1)  and  (4.2) 


PA  *  t 


22  21  b 


B 


'»  ‘  <rA-*U)  4  *"d 


■21 


(rB“i22)  bA 


(4.25) 


(4.26) 


b. 

b„ 

812l  "  lrA_8lll 

A 

bB 

and  | s21 | 

trB_822l 

B 

bA 

Having  found  s..  and  >22  from  equation  (4.3),  we  have  at  least  three  sets 
of  rA,  rB  and  of  |bA/bg|  using  (4.23),  since  no  ambiguity  of  sign  arises 
from  this  (because  the  sign  of  a  modulus  is  positive).  Hence  | s^2 1  aiu* 
|s2il  can  each  be  calculated  several  times  from  equation  (4.26)  and  their 
average  values  taken  (as  for  calculating  (KA/Kg)).  We  can  thus  measure 
| sx2 !  and  1 321 1  °£  a  two-port  junction  transmitting  the  same  mode  at  both 
ports  (whether  S12  ■  *21  or  812  *  s21^* 

4.5  DIRECT  MEASUREMENT  OF  ARGUMENTS  Lsl2  and  £s2i 

In  principle,  the  phase  angles  ^S]^  and  1.S21  can  be  measured  but  this 
process  has  not  been  implemented  at  RSRE  because  this  measurement  has  not 
been  required  and  because  the  network  analyser  of  Figure  3  has  to  be  modi¬ 
fied  to  include  isolators  (which  are  temperature-dependent  components  and 
which  would  limit  the  bandwidth  of  some  ref lectometers) .  Nevertheless,  we 
include  a  derivation  of  the  method, following  (19),  both  for  the  sake  of 
completeness  and  to  clarify  the  additional  assumptions  necessary  for  this 
measurement  to  be  possible. 

The  modified  network  analyser  would  take  the  form  illustrated  in 
Figure  6,  for  which  it  must  be  assumed  that  the  phase-shifter  is  precisely 


O. 


Source 


uV 
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& 

divider 
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•(see  text) 


FIG-6.  DUAL  NETWORK  ANALYSER  SUITABLE  FOR  MEASURING  ZS,?  AND 

^S,«  OF  PUT 
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repeatable  for  each  aeries  of  power  ratio  observations.  It  is  also  assumed 
that,  after  the  ref lactometers  have  been  calibrated  by  the  methods  already 
described  and  after  measurements  have  been  made  of  the  scattering  coeffi¬ 
cients  of  three  reciprocal  two-port  junctions,  then  the  flexible  waveguide 
is  replaced  by  a  rigid  one.  These  three  reciprocal  junctions  must  each  be 
of  the  same  length  £  as  the  non-reciprocal  two-port  DUT;  the  Z.«2l) 

of  two  of  them  must  differ  and  the  third  must  be  a f non-reflecting  line. 

This  line  must  be  capable  of  being  connected  at  XX  in  Figure  6  (ie 
between  the  power  divider  and  one  isolator)  so  that  the  connectors  forming 
this  reference  plane  must  be  of  the  same  form  as  those  of  the  ref lectometer 
measurement  ports. 

We  approach  the  measurement  problem  by  considering  the  entire  network 
analyser  as  a  three-port  junction,  of  which  two  ports  are  the  ref lectometer 
measurement  ports  and  the  third  is  a  reference  plane  P  situated  between  the 
power  divider  and  the  source,  as  shown  in  Figure  6.  We  can  write  the 
scattering  equations  for  this  junction  as 

8ab  8ap"  V 

8BB  8BP  aB 

*PB  ®PP-  -aP- 

and  the  first  two  of  these  equations  are 

bA  '  *****  *  ■*»*»  *  *AP*P  *”“  bB  *  *****  +  V>  *  *BP*P 
so  that  ap  can  be  eliminated  between  them  giving 

bA  •  *****  *  *AB*B  *  **Fa/*BP)(bB  -  “BA**  '  *BB*B> 

But  the  VRC  measured  by  the  ref lectometers  are  rA  ■  aA/bA  and  Tg  *  ag/bg 
so  that 


(A. 27) 


Mow  the  quantities  in  braces  in  equation  (A. 27)  are  characteristics  of 
the  network  analyser  and  if  they  can  be  found  by  calibration  then  bA/bg  can 
be  calculated  from  rA  and  Tg  for  each  phase-shifter  setting  allowing  *n 
and  S2i  to  be  completely  determined  from  equation  (A. 25). 

Calibration  of  the  network  analyser  for  measuring  l*\2  *nd  A*21  thu> 
involves  calibrating  the  repeatable  phase  shifter  and  the  rest  of  the 
analyser  and  it  is  necessary  for  this  that  any  change  of  phase  between  bA 
and  b«  caused  by  physically  separating  the  two  ref lectometers  be  known. 
Assuming  sufficient  isolation  by  the  two  isolators  of  Figure  6  for  the 
approximation  sAg  ■  sgx  ■  0  to  be  valid,  then  sm  and  sgg  will  not  change 
when  the  reflectionless  line  is  inserted  at  XX'  of  Figure  6  to  allow 
separation  of  the  ref lectometers  and  s^p  changes  by  the  factor  e~5l,  where 
6  la  the  known  propagation  constant  of  the  line  (as  does  SpA).  (This  can 
be  shown  by  analysing  the  network  of  power  divider,  isolators  etc  forming 
the  three-port  junction.) 
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Given  Chis  assumption  of  sufficient  isolation  and  the  approximations 
resulting  from  it,  the  two  ref lectometer  measurement  ports  can  be  connected 
to  each  other,  so  that  bA/bg  *rB,  and  we  have  from  (4.27)  that 


al  ♦  a2 


"  A3rB 


where 


{SAP/SBP} 


{SAA> 


tsBB8AP/sBPJ 


By  connecting  in  turn  the  two  known  two-port  reciprocal  junctions  between 
the  ref lectometer  ports  (for  which  we  also  need  to  connect  the  known  line 
at  XX'),  we  can  provide  from  equation  (4.25)  two  more  known  (bA/bg)m, 

(m  •  1,2)  so  that  from  (4.27)  we  have 

W.  -  *1  *  A2<VbB>„  rA»-A3  «"5tra. 

Hence,  we  can  obtain  a  set  of  simultaneous  equations  in  the  unknowns  Aj, 

A2,  A3  for  each  phase-shifter  setting.  Assuming  this  repeats,  this  further 
stage  of  calibration  allows  equation  (4.25)  to  be  used  to  obtain  S12  and 
S21  completely  for  a  non-reciprocal  two-port  junction. 

We  remark  that  values  of  s^  and  S21  so  obtained  are  most  suitable  for 
determining  Z.S32  and  L82l‘  The  magnitudes,  jsj^l  and  |s2]J» obtained  from 
(4.26)  are  likely  to  be  more  accurate  because  they  do  not  rely  on  the  phase- 
shifter  being  repeatable  or  on  the  isolation  being  sufficient.  Consistency 
between  the  magnitudes  obtained  by  (4.25)  and  (4.26)  may  provide  a  verifi¬ 
cation  of  these  assumptions.  We  also  remark  that  this  method  provides  a 
method  of  verifying  the  plausibility  of  arguments  used  in  section  4.3  for 
resolving  the  ambiguity  of  sign  of  the  square  root  of  (s^)2  for  a  recipro¬ 
cal  two-port;  even  with  the  perhaps  dubious  assumptions  made  concerning  the 
repeatability  and  predictability  of  the  three-port  junction  comprising  the 
analyser  of  Figure  6,  the  method  should  be  able  to  resolve  an  ambiguity  of 
0  or  » . 

4.6  SIMULTANEOUS  COMPARISON  OF  POWER  STANDARDS 

We  have  shown  in  section  3.4  that  power  standards  can  be  compared  by 
connecting  them  in  turn  to  the  measurement  port  of  a  six-port  ref lectometer 
and  using  equation  (3.22).  We  note  that  a  dual  ref lectometer  network 
analyser  allows  equation  (4.24)  to  be  used  for  comparing  two  power  standards 
simultaneously  (rather  than  sequentially)  and  their  VRCs  can  be  also 
determined. 

5  AIR  DIELECTRIC  TRANSMISSION  LINES  AS  IMPEDANCE  STANDARDS 

The  assumptions  made  in  our  theory  for  calibrating  the  dual  reflectometer 
are  that: 


a.  the  source  provides  single-frequency  radiation; 

b.  each  detector  provides  an  indication  proportional  to  the  power  it 
absorbs ; 

c.  the  reflectometer  junction  transmits  only  one  mode  at  each  port  and 
that  both  measurement  ports  transmit  the  same  mode; 
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d.  the  uniform  line  forming  the  calibration  standard  has  measurable 
attenuation  and  connectors  at  each  end  of  a  type  that  permit  its  length 
to  be  accurately  measured. 


Then,  to  the  extent  that  the  connectors  at  the  ends  of  the  standard  line  may  be 
considered  electrically  "invisible*1,  both  ref lectometers  (including  the  connec¬ 
tors  forming  their  measurement  ports)  are  calibrated  in  terms  of  the  impedance 
Zc  characterising  the  mode  transmitted  by  the  standard  line.  Connectors  to 
realise  condition  (d)  for  rectangular  waveguide  are  relatively  easily  implemented 
by  lapped,  dowelled,  waveguide  flanges;  and,  for  some  sizes  of  air-dielectric 
precision  coaxial  line,  connectors  are  available  providing  a  good  approximation 
to  coplanar  inner-  and  outer-conductor  contacts.  For  both,  the  contact  impedance 
is  included  in  the  attenuation  constant  6  measured  of  the  standard  during  cali¬ 
bration,  but  this  inclusion  is  unfortunately  not  self-consistent  (because  when 
the  measurement  ports  are  connected  together  only  one  connector  pair  is  involved, 
whereas  connection  through  the  standard  line  involves  two  pairs). 

At  some  frequencies,  the  imperfect  surface  finish  of  coaxial  line  conductors 
causes  Z0  to  depart  measurably  from  the  nominal  characteristic  impedance  R0  (real) 
calculable  for  the  TEM  mode  for  an  ideal  lossless  line  of  the  same  dimensions. 

In  this  event,  the  measured  VRC  Tx  and  Fg  may  be  renormalised  to  the  nominal 
characteristic  impedance  R0.  The  transmission-line  model  of  a  two-conductor 
line  describes  it  in  terms  of  series  resistance  R  and  inductance  L,  shunt 
conductance  G  and  capacitance  C,  all  per  unit  length  at  some  particular  angular 
frequency  u;  then 

Zq  ■  /  (R+ju>L)  /  (G+ jwC)  and  6  ■  /(R+jtuL)  (G+jwC) 

where  6  may  be  identified  as  the  true  propagation  constant  of  the  line  as  defined 
in  section  3.3  and  ZQ  is  the  impedance  characterising  the  line.  The  branches  of 
the  square  roots  chosen  are,  of  course,  those  which  make  ZD  positive-real  and  6 
positive-imaginary  for  small  R  and  G  so  that  x  ■  exp (-261)  in  section  3.3.  Hence 
on  rearranging  we  have  ZQ  ■  5/ (G+juC)  where  C  can  be  identified  as  the  calculable 
electrostatic  capacitance  per  unit  length  of  the  precision  line.  For  air 
dielectric  we  may  accurately  set  G  ■  0  so  that  ZQ  may  be  determined  from  the 
measured  6  and  the  calculable  C.  Moreover 


where  vQ  is  the  velocity  of  propagation  so  that  the  indicated  T  (normalised  to 
Z0)  may  be  renormalised  to  Rc  (to  give  TQ)  since 


l  ±2- 

o  l-r 


l+r 


R 


o  i-r 


ie 


r (l+Z  /R  )- (1-Z  /R  ) 
-  m  o  o  oo 

°  (l+Z  /R  )-r(l-Z  /R  ) 
o  o  o  o 


Although  this  process  of  renormalisation  to  a  standard  R.  calculable  from 
mechanical  dimensions  and  the  velocity  of  propagation  in  air  is  attractive  as  a 
means  of  standardisation,  it  will  require  subsidiary  measurement  of  the  "residual" 
contact  impedance  of  the  air-line  connectors  as  this  enters  into  the  propagation 
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constant  measured.  The  resolution  o£  the  six-port  ref lec tome ter s  should  enable 
these  to  be  applied  in  established  methods  of  measurement  of  this  contact 
impedance  (16) . 

6  SIMULTANEOUS-EQUATION-SOLVING  ALGORITHMS 


6.1  LEAST-SQUARES  SOLUTION  FOR  REAL  UNKNOWNS 


This  section  outlines  the  method  employed  to  refine  the  provisional 
solutions  of  simultaneous  non-linear  equations  appearing  in  the  six-  to 
four-port  reduction  derived  in  section  3.1.  The  methods  are  outlined  here 
because  most  desktop  computers  have  no  library  of  suitable  routines.  One 
stage  of  the  method  entails  solving  an  overdetermined  set  of  linear  equa¬ 
tions.  The  solution  of  these  is  discussed  separately  here  because  the  same 
problem  occurs  in  calibrating  either  a  four-  or  a  six-port  ref lectometer  by 
the  use  of  known  loads,  derived  in  sections  2.2  and  3.2,  respectively. 


Suppose  we  have  a  set  of  M  equations,  not  necessarily  linear,  in  N 
(<  M)  unknowns;  thus 

fi(x1,  x2,  . ...  Xj,  ...  Xjj)  -  0  for  i  -  1,  ...»  M  (6.1) 

If  we  have  approximations  °x.  to  the  unknowns,  we  may  re-write  (6.1), 
neglecting  second-order  terms,  as 


N 


f i (°xi ,  °x2  ....  °xN)  +  ^  |  (JXj  -  °x.)  ■  0  for  i  ■  1,  ...,  M  (6.2) 

j-1  j 

where  all  the  partial  derivatives  are  evaluated  at  the  point  (xj  ■  °x^,  .. 
xj  ■  °xj.  ...,  xjj  -  °xm).  This  set  of  equations  is  linear  in  the  correc¬ 
tions  (ixj  -  °Xj) ,  and  may  be  solved  for  them  by  the  method  described 
below.  The  correction  process  may  be  repeated  as  often  as  necessary; 
it  displays  quadratic  convergence,  like  Newton's  method  of  which  it  is  a 
generalisation. 


Suppose  now  that  we  have  a  set  of  M  linear  equations  in  N  (<  M) 
unknowns;  in  matrix  notation  we  may  write 

Ax  -  b  (6.3) 


where  A  contains  the  coefficients  and  is  of  order  (M  x  N),  x  contains  the 
unknowns  and  is  of  order  (N  x  1) ,  and  b  contains  the  right-hand  sides  and 
is  of  order  (M  x  1).  We  assume  for  the  moment  that  A,  x  and  b  are  real; 
this  restriction  will  be  removed  later. 


In  the  special  case  where  M  ■  N,  (6.3)  admits  the  solution 

x  -  A_1b  (6.4) 

When  M  >  N,  it  is  customary  to  choose  the  solution  x  such  that  the  scalar 
quantity  (Ax-b)T(Ax-b)  is  an  absolute  minimum,  where  T  denotes  matrix 
transpose.  This  is  the  "method  of  least-squares",  first  adopted  by  Gauss 
in  the  light  of  his  error  distribution  but  now  routinely  used  because  of 
its  formal  simplicity. 
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To  minimise  (Ax-b)T(Ax-b)  absolutely,  we  differentiate  partially  with 
respect  to  each  component  of  x  and  set  all  the  resulting  quantities  equal 
to  sero.  We  then  obtain  N  simultaneous  linear  equations  in  the  N  components 
of  x,  and  they  admit  the  formal  solution 

x  -  (ATA)_1ATb  (6.5) 

When  M  •  N,  (6.5)  reduces  to  (6.4),  the  absolute  minimum  of  (Ax-b)^(Ax-b) 
then  being  sero. 

When  M  >  N,  (6.5)  is  not  numerically-stable;  the  three  appearances  of 
A  (compared  to  the  one  appearance  in  (6.4))  cause  a  large  increase  in  the 
rounding  errors  during  the  actual  solution.  This  problem  is  overcome  by 
resolving  the  matrix  A  into  a  product  of  matrices  Q  and  U,  where  the  matrix 
Q  (of  order  (M  x  N))  has  the  property  that  QTQ  -  %  (the  identity  matrix  of 
order  (N  x  N)),  and  the  matrix  U  is  of  order  (N  x  N)  and  can  be  inverted 
easily.  From  (6.5)  we  then  have 

x  »  U_1QTb  (6.6) 

which  is  a  numerically-stable  form  for  x  because  the  components  of  A  only 
appear  once  instead  of  three  times  as  in  (6.5).  The  resolution  is  brought 
about  by  applying  to  A  N  successive  transformations  of  the  kind  named  after 
Householder,  which  are  described  by  Fox  (17).  The  matrix  U  is  upper- 
triangular  (it  has  non-zero  elements  only  on  and  above  the  diagonal  from 
its  top  left  corner  to  its  bottom  right  corner)  which  makes  it  easy  to  find 
x  from  (6.6)  or  from  the  equivalent  form 

Ux  -  QTb  .  (6.7) 

During  the  resolution  it  is  usually  necessary  to  alter  the  order  of  the  rows 

of  A;  this  process,  and  the  later  unscrambling  of  it,  will  also  be  found  in 

Fox' 8  book. 

6.2  EXTENSION  TO  COMPLEX  UNKNOWNS 

We  consider  now  the  extension  of  this  method  to  complex  A,  x  and  b. 

If  j  is  /-T,  we  may  write  in  an  obvious  notation 

A  -  Ar  +  jA.  ;  x  -  xr  +  jxj^  ;  b  -  by  +  jb.  (6.8) 

Then,  from  (6.3)  and  (6.8),  we  have 

(Arxf  -  AjX^)  ♦  j(Arx.  +  Aixr>  "  <br>  +  i<bi>  (6.9) 

If  we  equate  real  and  imaginary  parts  in  (6.9),  we  may  re-write  it  in  the 
partitioned -matrix  form 


which  is  purely  real  and  may  be  solved  using  the  method  of  the  previous 
section.  In  the  real  case,  minimising  (Ax-b)^(Ax-b)  is  the  same  as 
minimising  the  sum  of  the  squares  of  the  differences  between  the  left- 
and  right-hand  sides  of  each  equation; 
when  a  complex  problem  is  attacked 
using  (6.10),  we  minimise  the  sum  of 
the  squares  of  the  real  and  imaginary 
parts  of  the  differences,  which  is  the 
sum  of  the  squares  of  their  moduli. 

For  uncorrelated  equations  this  can  be 
shown  to  be  correct  from  probability 
considerations.  It  is,  apparently, 
unsatisfactory  because  the  phases  of 
the  differences  are  not  constrained  by 
minimising  the  sum  of  squared  moduli; 
but  by  considering  Figure  7  we  can  see 
that  this  does  not  matter.  We  have  a 
right-hand  side  vector  b  added  to  a 
difference  vector  5;  if  | 6 |  <  |b|, 
then,  wherever  the  direction  of  6  may 
be,  the  sum  vector  s  is  still  tightly-constrained  both  in  magnitude  and  in 
direction.  Moreover,  the  maximum  difference  in  phase  between  s  and  b 
corresponds  to  the  minimum  difference  in  magnitude,  and  vice  versa,  so  the 
differences  are  well-balanced  in  practice. 

7  CONCLUSION 

It  is  axiomatic  that  the  quality  of  a  product  or  service  cannot  be  main¬ 
tained,  let  alone  improved,  without  means  of  measuring  that  quality.  This  is 
true  also  for  measurement  processes  which  require  repeated  (and  therefore  in 
practice  automated)  measurements  for  their  variation,  and  the  stability  of  the 
reference  standards,  to  be  quantified.  In  principle,  the  calibration  methods 
treated  in  this  report  allow  the  stability  of  six-port  ref lectometer  instrumen¬ 
tation  to  be  checked  without  reference  to  standards  whose  long-term  stability 
would  have  to  be  inferred  or  even  assumed.  Moreover,  the  resolution  obtained 
at  RSRE  of  measurement  of  VRC  with  these  methods  suggests  that  the  long-term 
stability  of  the  reference  standards  can  be  compared  since  each  calibration 
requires  only  one  such  standard. 

The  derivation  of  the  theory  of  these  six-port  ref lectometer  calibration 
and  measurement  methods  has  been  presented  in  this  report  as  an  analytic  deduc¬ 
tion  from  the  basic  postulates  of  linearity  of  the  ref lectometer  and  its 
detectors.  This  was  a  necessary  prerequisite  for  programming  the  solutions  to 
the  problems  (which  are  trivial  only  in  the  mathematical  sense  that  they  have 
been  solved),  especially  to  ensure  that  all  the  ambiguities  of  square-root  sign 
are  resolved.  This  approach  appears  to  the  authors  to  be  more  appropriate  for 
computation  than  that  contained  in  earlier  treatments  (relying  on  geometric 
interpretations  of  diagrams  in  the  complex  plane  of  six-port  ref lectometer 
behaviour) . 
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